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Possible signatures of a super- Tonks-Girardeau gas in bosonic systems of trapped quasi-one- 
. dimensional dipoles are discussed at zero temperature. We provide estimation of the frequency of 

' the lowest compressional mode and compare it to analytical results derived using harmonic approach 

, in the high density regime. We construct an exact mapping of the ground-state wave function of 

. one-dimensional dipolar system of bosons, fermions and Bose-Fermi mixture and conclude that local 

properties and energy are the same at zero temperature. A question of to which extent the dipolar 
^j;^ I potential can be treated long- or short- range is discussed. 

oo ■ 

CN ; I. INTRODUCTION 

JrJ , Recent realization of chromium condensates pj] have brought much interest to unusual properties of dipolar con- 
(— I ' densates. Use of chromium have proved to be very advantageous as it has very large permanent dipolar moment of 
> ' /i ~ 6fiB- Other possible realizations of dipolar gases includ polar molecules [z, 's', '3| and laser-induced dipoles^, '6|. 
^ ■ Study of collective oscillations is a very sensitive tool for the investigation of the cold gas properties [7j]. Dipolar forces 
compete with short-range s-wave scattering interaction. Although in the first experiments the s-wave interaction was 
giving dominant contribution to the energy, the effects of dipolar interactions have been detected in the shape of 
■ expanding cloud [8]. During the expansion of an initially trapped gas dipolar forces might lead to an anisotropic shape 
. of an expanding cloudQ. Use of Feshbach resonance technique have brought experiment to a new level. Indeed, by 
' the means of Feshbach resonance technique it is possible to reduce s-wave isotropic contact interaction, such that the 
O I anisotropic magnetic dipole-dipole interaction between ^^Cr atoms becomes comparable in strength 10]. This induces 
, ^ , ■ large changes of the aspect ratio of the cloud, and, for strong dipolar interaction, the inversion of ellipticity during 
expansion can even be suppressed [Toj. Chromium atoms have have as well fermionic ^■^Cr isotope and Bose - Fermi 
cold mixtures have been recently realized [ll|- 
^ ' Furthermore, dipolar interactions are very interesting as they contain both repulsive and attractive parts. Having 
attractive interactions between bosonic particles at temperatures so low that condensate is formed might lead to 
J— ^ ' instability and a violent collapse process ITsl . IT^ . Collapse of large dipolar condensates in traps have been studied 
using non-local Gross-Pitaevskiifisl [l6l. llSj theorv and more precise Diffusion Monte Carlo (DMC) [13] approach. 



o 
o 



X 



One-dimensional (ID) cold systems have received great attention in the last years |T1, [S^, [2lj. Role of quantum 
fluctuations is increased in reduced dimensionality leading to sometimes very different behavior. For example, a 
peculiarity of a one-dimensional world is an absence of a true Bose condensate in ID homogeneous systems even at 
zero temperature (22j . There is a certain trend to study low-dimensional systems in the last years. 

One of the most precise techniques, which can be used for testing the equation of state, is the measurement of 
the frequency of the "breathing" mode produced by a sudden change of the frequency of the harmonic trapping. We 
note that differences of several percent can be resolved in present high-precision experiments with cold gases (see, for 
^ • example, Ref. 23]). 

5^ ] In this paper we address properties of a quasi-one-diniensional dipolar system at zero temperature. We provide an 
explicit mapping of the fermionic ground-state wave function and a wave function of Bose-Fcrmi mixture to bosonic 
wave function in one-dimension. Due to this mapping energy of dipolar systems containing fermions is predicted using 
results previously obtained for bosons. We calculate frequency of the lowest breathing mode in the trap and compare it 
to analytic predictions obtained in the high density limit. Measurement of the frequencies of collective oscillations can 
provide a signature of a super- Tonks-Girardeau (STG) regimefU, [2^ described by very strong correlations (stronger 
than in the Tonks-Girardeau gas (26j in which the coupling constant is infinitely large). It is very difficult to reach 
STG regime in systems with short-range interactions, while its realization is feasible in dipolar systems. We address 
problem of superfluidity in one-dimensional systems and a question up to which extent dipolar interaction is short- 
er long- range. 
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II. MODEL 

A system of N dipoles in one-dimensional geometry is described by the following model Hamiltonian: 

i—l * i—1 i<j ' ' 

Here we assume that all dipoles are polarized and are oriented perpendicularly to the one-dimensional line. This 
stabilizes the system (as only repulsive part of dipolar interactions is relevant) and avoids collapses due to attraction. 
Expression for the coupling constant Cdd depends on the nature of the dipolar interaction. For example, following 
realizations are possible: 

1) Cold bosonic atoms, with induced or static dipole momenta, in a transverse trap so tight that excitations of the 
levels of the transverse confinement are not possible and the system is dynamically one-dimensional. The longitudinal 
confinement is described by the frequency uj^ of the harmonic trapping potential. The dipoles themselves can be either 
induced or permanent. In the case of dipoles induced by an electric field E the coupling has the form Cdd — E'^o? ^ 
where a is the static polarizability. For permanent magnetic dipoles aligned by an external magnetic field one has 
Cdd = '71^, where m is the magnetic dipole moment. We suppose that s-wave scattering length is tuned to zero by 
applying Feshbach resonance and only dipolar forces are relevant. 

2) Spatially indirect excitons in two coupled quantum wires. A quantum wire is a semiconductor nanostructure 
where an electron or a hole is allowed to move only in one direction and excitations of the transverse quantization 
levels are negligible. In two parallel quantum wires, one containing only holes, and the other only electrons, holes 
and electrons couple forming indirect excitons. If such a system is dilute enough, it constitutes a \D set of dipoles. 
In this case Cdd = e^d? /e, where e is an electron's charge, e is the dielectric constant of the semiconductor, and d is 
the distance between the centers of the quantum wires. This system is ID counterpart of 2D indirect exciton system 
in coupled quantum wells. 

III. SUPER-TONKS-GIRARDEAU REGIME 

Properties of a homogeneous system (wz = 0) have been studied numerically infill by means of Diffusion Monte 
Carlo method. The equation of state and correlation functions have been calculated as a function of the guiding 
parameter nro, where — MCddl {^i^'h^) is a characteristic length, n being the linear density. It has been found 
that the system is extremely correlated and shows crystal-like properties as the strength of dipolar interactions nrg is 
increased. This behavior is very different from the one of systems with short-range potentials, like the ones recently 
realized experimentally [1^ [2^. Indeed, one-dimensional systems with (5-interaction, Viz) = —2U^/maiDS{z), (Lieb- 
Liniger model[30|]) have a completely different behavior of correlation functionspBll. |33 | which in this case is much 
closer to the one of a weakly interacting Bose gas. For example a typical shape of the static structure factor S{k) 
is a smooth function which goes from zero for fc = to asymptotic constant value S{k) = 1 for large In the 
most strongly interacting limit of Licb-Liniger system —2h^/maiD +oo (Tonks-Girardeau regime |26j) the static 
structure factor is equal to the one of an ideal Fermi gas and has a discontinuity in the derivative at momentum 
fc| = 2fci? = 27rn. Even stronger correlations might be achieved in a system with short-ranged potential by quickly 
crossing the confined induced resonance (33.] as proposed in (23 |. In this resulting "super- Tonks-Girardeu" regime 
bosonic atoms will interact with attractive potential —2T%^/mai]j —oo. The true ground state in this regime is 
a soliton-like state with large and negative energy [3l| ("attractive collapse" of the system). Instead, if the sweep 
across resonance is fast and does not significantly perturb position of the particles, the system will still remain in 
a gas-like state (which is metastable). As shown in Ref. '24!|, this state is dynamically stable if the gas parameter 
naiD is relatively small. The static structure factor in super- Tonks-Girardeau regime has a peak at \k\ = 2kp and 
height of the peak increases for larger values of the gas parameter nan). The super-Tonks-Girardeau regime also 
exists in a number one-dimensional system (bosonic and fermionic) with infinite strength potentials, namely hard- 
rods [1^ and Calogero-Sutherland model 36]. So far this regime have never been observed experimentally with the 
best candidate for its observation being a system of dipoles. An experimental signature of STG regime in a trapped 
system is a frequency of lowest compressional mode ilz larger than 2lUz [24] . In the Section IIVI we calculate explicitly 
the dependence of Qz on parameters of a trap. 

Super-Tonks-Girardeau regime is expected to have quite unusual properties. It has been shownjs^, IH, [3§] for a 
Lieb-Liniger model that the dynamic form factor has a power low in the point where the excitation spectrum touches 
zero S{lu, 2kp) oc w''^^, 0, where ry — 2hkF/'mc. In the regime of repulsive J-interaction r] > 2 and 5(0, 2fci?) = 0. 
For the marginal case of Tonks-Girardeau regime 5'(u;, 2kp) = const. Instead for attractive (5-interaction (i.e. in the 
Super-Tonks-Girardeau regime) this expression predicts a weak (power law) divergence in the dynamic form factor 
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FIG. 1: Square of the lowest breathing mode frequency, £7^, as a function of the coupling strength Nr^/a^^. The solid line: using 
equation of state of one-dimensional dipols (data is taken from |49l] . Fig. 6.5), dashed line: high density expansion, Eq. (|3]). 



5'(w, 2kp) 00, w — >■ 0. For additional information on singularities in the dynamic structure factor for different values 
of momenta see Ref. [13, EH 

The problem of superfluidity has peculiarities in one-dimensional system. While different ways to define superfluid 
part in three- and two- dimensional systems are consistent, this is not the case in a one-dimensional world. Indeed, its 
calculation as a response of a liquid to sample boundary motion [winding-number method [42]) for the highly correlated 
states described by translationary invariant gas-like wave function (for example, exact wave function for the Tonks- 
Girardeaujl^, Calogero-Sutherlandfisl l43| . hard-rods systems; and DMC evaluation for Licb-Liniger[3l|, dipolar 
[131 systems) would provide the result that such systems are completely superfluid. This argument would apply even to 
the TG regime, where bosonic system has many similar properties to an ideal Fermi gas. At the same time, one should 
keep in mind that exposure of an one-dimensional ideal Fermi gas to a tiny perturbation will change the ground-state 
wave function in a dramatic way: the overlap of the new wave function and the old is essentially zero. This effect 
known as orthogonality catastrophe (see, for example, textbook ^45'|) shows that in one-dimensional system stability to 
external perturbations have to be carefully checked. Contrary to winding number approach. Landau argument would 
lead to exactly opposite result classifying systems as completely normal. Indeed, in one-dimensional Luttinger liquids 
the excitation spectrum always touches zero at finite value of momentum, \k\ = 2kp (see, for example, Ref. [39]; for 
Lieb-Liniger model the "type 11" excitation[46] that touches zero can be identified as a dark soliton^^)- Another 
way to calculate superfluid density by the response to transverse probe (transverse current-current response) is not 
applicable within one-dimensional description as no transverse direction is included in the model. Contradictions 
in the results arise from different definitions of superfluid part and reflect the non-standard nature of the system. 
Probably the most natural and appropriate way to test the superfluidity in ID world is done by dragging a small 
impurity (perturbation) through the system and seeing if this leads to energy dissipation. The force Fy, experienced 
by the system, depends on the interaction parameter rj as Fy oc y'^^[38] where V is velocity with which a small 
(5-perturbation moves through the system. Thus in the mean-fleld limit ?7 — > c» the force is vanishing and from a 
practical point of view system behaves analogously to a superfluid. On the contrary, in the Tonks- Girardeau limit 
Fv oc V and system behaves, from the point of view of friction, as a normal system, where the drag force is proportional 
to the velocity. In between there is a smooth crossover. The dipolar one-dimensional systems are expected to behave 
as normal ones. 

Although the Bose-Einstein condensation is absent in a one-dimensional system even at zero temperature, its 
reminiscence still can be observed in divergence of the momentum distribution for fc = 0. This divergence is present 
in a homogeneous system of bosons in the Tonks-Girardeau regime. From the Luttinger liquid theory it is possible to 
show that deep in super- Tonks-Girardeau regime this divergence will be removed. This happens for rj = I (similarly 
to Calogero-Sutherland [1^ and Hard-Rod [33 systems) or in terms of Luttinger parameter K — 1/2. The Luttinger 
parameter can be extracted from the equation of state or from the phononic part of the static structure factor and is 
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IV. FREQUENCIES OF COLLECTIVE OSCILLATIONS 

The equation of state obtained in [27,] can be used to predict properties of a trapped one-dimensional gas within 
local density approximation (LDA) (see, for example, [50| ) by assuming that the chemical potential in a trap is a sum 
of local chemical potential, taken to be the same as in a homogeneous system fJ-homin), and the external harmonic 
potential: ^ — fJ.hom{n{z)) + (l/2)Ma;^z^. Within LDA number of particles N and oscillator length az = y/h/Muj^ 
come in a single combination, Nr^/al. The LDA is expected to describe correctly properties of a gas in the trap if 
size of the cloud is large compared to az. From the knowledge of the density profile n{z) one can obtain the mean 

square radius of the cloud (z^) = J^j^n{z)z'^dz/N and thus, making use of the result [sJI fi^ — —2{z'^)/{d{z'^)/duj1) 
to calculate the frequency il^ of the lowest breathing mode. 

Fig. [1] shows the square of the lowest breathing mode frequency Sl^ as a function of the coupling strength Nr^jd}^. 
In the TG regime, Nr\ /a1 <^ 1, the frequency reaches a typical result of an ideal Fermi gas Q,z — 2lUz. The observation 
of a breathing mode with a frequency larger than 2wz would be a clear signature of the super- Tonks-Girardeau regime. 

Attempt to calculate analytically the equation of state in the regime nrp <C 1 in terms of the s-wave scattering 
length and to obtain expansion for the frequencies of collective oscillations encounters problems as will be explained 
in Section IVll Instead some analytical results can be obtained in the regime of high densities nro ^ 1. In this regime 
the dipolar system has crystal-like properties [13]. At the same time the difference in the energy between gas-like and 
crystal-like descriptions is extremely small. This justifies an attempt to derive in a perturbative way the equation of 
state of a crystal and use it for an approximate description of the gas-like phase at the same density. We use classical 
crystal harmonic approach to calculate equation of state in the perturbative manner. Details of the calculation are 
provided in Appendix [X] The result for the equation of state is 



S - C(3)^Ko)3 + C^Ko)^/^ + ...,nro » 1, (2) 



with C ~ 2.26.... Application of the LDA to "perturbative" equations of state has been studied in details in [50|. 
Density profile, total and release energy are easily obtained by using expansion ^ in the high density Nr^/a^ ^ 1 
regime. Here we report the frequency of the lowest compressional mode: 

1^ ^ 1.1358... 

The obtained dependence is shown in Fig. [1] as a dashed line and provides good description for large densities. 



V. BOSE-FERMI MAPPING FOR THE GROUND-STATE WAVE FUNCTION 

In the previous paper [13] we have calculated zero-temperature equation of state for the homogeneous Hamilto- 
nian ([T]) with = for bosonic particles and have sampled correlation functions over the ground-state wave function 

b{zi, zn)- As in experiments dipoles can be not only bosons, but as well fermions, or even a mixture of bosons 
and fermions we note that the equation of state we have obtained is applicable also to systems containing fermions. 
In order to prove that we construct an exact mapping of the wave function of bosonic dipoles to a wave function of 
fermionic wave function by analogy to what Marvin Girardeau did in his classical work[26j. For simplicity we start 
with a system of same-spin fermions. Such a system is described by the Hamiltonian ([1]). The main difference from 
the bosonic case is that fermionic wave function must be antisymmetric with the respect of exchange of any two 
particles. This can be done as 

^f{zi,-.-,zn) = ]^sign(z^ - Zj)^B{zij ■■■,zn)- (4) 

It is easy to check that the symmetry of wave function ([4]) is correct. Furthermore, due to Pauli exclusion principle, 
two Fermions are not permitted to stay in the same place f{zi, zat) = if = Zj. This is already satisfied in the 
construction of bosonic wave function ^ b{zi, zpf) due to divergence of l/jz]^ interaction for small z. It means that 
performed Diffusion Monte Carlo calculation for bosons is equivalent to Fixed-Node Diffusion Monte Carlo calculation 
for fermions with exactly known nodal structure as far as energy and local quantities (pair correlation function, static 
structure factorjete.) are concerned. The trick Q has been successfully used for study Fermionic Calogero-Sutherland 
model in Ref. [39]. 
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Another system which will have similar zero-temperature equation of state is a system of two-component fermions, 
where the mass of different spin atoms a =1, J, is the same. Hamiltonian of such a system is given by 

H=-j^ E ^+^^^^-n:-^.^+^i>:^r— r^+>:^-^+>:^r— ^1 (5) 

Ground state wave function of system with totally N — iV^ -I- N-^ atoms is than mapped onto a system of N bosons 
with a Hamiltonian ([T]) as 

^f{zI, ...,z'^j^,,zi,...,zl^i) = Y[sign{zj - zj) ]Jsign(z| - z^^)^ b{zI, z'^^, z^, z^j^). (6) 

i<j i<3 

Thus we conclude that a system one-dimensional Fermionic dipoles has the same ground-state equation of state 
as a system of bosonic dipoles. Consequently the LDA density profile in a trap is the same. This means that the 
frequencies of the lowest breathing mode for fermions follow the dependence shown in Fig.[l]with N being total number 
of dipoles. The limit Nr^/a?^ — > corresponds to ideal fermions and ideal fermions has the spherical breathing mode 
il — 2u}fio in any dimension (see, for example, [l^l). For a finite value of Nr^/al the frequency is increased due 
to repulsive interactions. Similar effects have been predicted in Ref. [5^ using Bethe-ansate theory for systems of 
two-component fermions with J-pseudopotential attractive [H, [s^l and repulsive [5^1 interaction between atoms of 
different spin. While homogeneous system of bosons collapses if the interaction is attractive, this is not the case for 
fermions, where Pauli priiiciple stabilizes the system. Adding weak attraction between atoms leads to softening the 
breathing mode fiz < 2ujz [5^ . 

The Bethe ansatz method permits to find ground-state energy and, thus, study collective oscillations within LDA in 
a mixture of one-dimensional bosons and fermions with i5-pseudopotential repulsive interactions and arbitrary bosons- 
fermions density ratiojsl]. Such a mixture was found to be always stable against demixingjs^. Ground-state wave 
function of a mixture of one-dimensional bosons and fermions can be obtained with the same reasoning as for Q 
and ([n]), thus leading to the same frequency of oscillations as in Fig. [T] Such kind of mapping turns out to be quite 
general. Indeed, it works for one-dimensional systems where one-dimensional interaction potential diverges when two 
particles meet and shows that energy, local quantities (pair-correlation function, three-particle correlation function, 
etc.) are the same for bosons and fermions at zero temperature. 



VI. ARE DIPOLAR INTERACTIONS LONG-RANGE? 



It is common to oppose long-range dipolar interactions to short-range interactions described by s-wave scattering 
length. But are dipolar interactions really long range? Or to which extend are they long-range? There are at least 
two ways to approach this question. 

The first way is to classify potential Vint (^) as long- or short- range depending if the chemical potential is extensive on 
intensive quantity. If the homogeneous properties of a large system can be defined by density n only (in appropriate 
units), fj, = fJ.{n), the potential is short-range. If instead number of particles N have to be explicitly specified, 
H = fi{N,n), due to strong (diverging) dependence on N, the potential is of a long-range. This can be immediately 
checked by testing the convergence of the potential energy at large distances: 

oo 

1= J V„u(r)r''-Ur, (7) 

where Lmin is some cut-off length and D is dimensionality. Following this definition potential is short-range if it 
decays at large distances faster than r"^ in D dimensions (see, for example, ^7}). From this point of view, l/|rp 
potential is long-range in 3D, while it is short-range in ID and 2D. 

Alternatively, short-range potentials can be defined as potentials that can be described by the asymptotic phase 
shift. This means that a short-range potential of range R can be approximated at large distances r ^ i? by a free- wave 
with an appropriate phase shift, or, being the same, the s-wave scattering length a. At sufficiently small densities 
na^ — > 0, the only relevant length is a and properties (for example energy) can be expressed in terms of the gas 
parameter na^ . We use definition of the s-wave scattering length as a position of the node of analytic continuation of 
the scattering solution from distances larger than the range of the potential in the zero-energy scattering limit. This 
definition works well in three-dimensional systems, but also it is applicable to low-dimensional systems. 
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It is possible to solve the two-body scattering problem for potential at zero energy and find the scattering 

solution /(r). One needs to look for a regular solution (/(O) — 0) of the following differential equation: 

d'fjr) D-ldfir) fir) 

a"l T ' 3~ ^ ^' voj 

or^ r or r-^ 

where length is expressed in units of tq and energy in /Mr^. Solutions of ^ can be written explicitly f^^{r) (x 
r^''^Ki{2r-^/^), f^{r) cx i^o(2r-i/2)^ /^^(r) cx r-^/^ Ki{2r-^/^), where Knir) denotes the modified Bessel function 
of the second kind. In order to find the s-wave scattering length one has to expand /(r) far from the range of the 
potential (r — > oo) and compare it to the similar expansion of a plane-wave in appropriate number of dimensions: 
ffreei'^) oc 1 - a^D / r , fj^eei^) °^ ln(r/a2D), ffreel^) oc T - fliD- Expanding the solutions of Eq. ^ we find 
f^{r) cx 1 + (27 - l)/r - (lnr)/r + ©(r'^), p°{r) cx 27 - Inr + 0{r-^), p^{r) cx r -h (27 - 1) - Inr -f 0{r-^), 
where 7 = 0.577... is Euler's constant. There are logarithmic terms appearing in all dimensions. In a two dimensional 
system such a term is compatible with asymptotic behavior of a free- wave solution and it is possible to define a finite 
scattering length 020 — e^^rp — 3.172...ro. Indeed, numerical evaluation of the equation of state in a two-dimensional 
dipolar systemfs^ is in agreement with equation of state for hard-disks with same values of a9n[59l|. Also, Bogoliubov 
theory at small rta|^ provides correct predictions for the correlation functions and condensate fraction. Thus, in 
two-dimensions dipolar l/lrp potential can be treated as a short-range one. Of course, values of the typical value 
of the gas parameter na^D for which universal in terms of na2D description starts being valid is smaller for dipolar 
interaction, than for "usual" short-range potentials (for example, a soft-disks[59|). 

Instead, in a one-dimensional system the free-wave solution does not contain logarithmic term and this makes the 
dipolar solution be incompatible with it. For example, for short-range potentials s-wave scattering length can be 
calculated through the limit lim [r — f{r)/f'(r)] ~ an). As it is easy to see, corresponding expression diverges in the 

r — >QC 

case of dipolar interaction potential. It means that it is not possible to describe properties of a one-dimensional dipolar 
system by a short range potential. The same conclusion can be reached by comparing the ground-state energy of a 
repulsive dipolar system [23, SI] to ground-state energy of the repulsive short-range 5-pseudpotential S^l in the dilute 



regime noiD <C 1. The leading term in the energy E is the one of the Tonks-Girardeau gas Eth/N = ir'^h'^ n\j^ I &M\2(^. 
The first correction should include (if the short-range description is possible) terms naiD- Instead it is clear that it is 
not the case, as for Lieb-Liniger gas aio < and energy is lowered E < Etg- By simply neglecting the logarithmic 
term for dipoles one would find negative scattering length am = (1 — 27)ro = —0.154..., and the energy would be 
lower. But we definitely know that this is not the case, instead energy is higher E > Etg- Thus we find problems 
trying to describe the dipolar interaction potential by a short-range model. Most probably, the same will happen in 
a three dimensional system. 



VII. CONCLUSIONS 



In conclusions, we have studied possible signatures of a super- Tonks-Girardeau gas in a system of trapped quasi-one- 
dimensional dipoles at zero temperature. This regime can be entered by exploiting a confinement induced resonance 
of the effective ID scattering amplitude. Using previously calculated equation of state we provide predictions for the 
frequency of the lowest compressional mode. Properties in the high density regime are calculated within harmonic 
approximation. We provide an explicit mapping of the ground-state wave function of one-dimensional dipolar system 
of bosons; fermions and Bose- Fermi mixture and conclude that local properties and energy are the same at zero 
temperature. A question to which extent the dipolar potential can be treated long- or short- range is discussed. 
Superfluidity in one-dimensional systems is tested using different definitions. 

During preparation of the paper to publication a related article [60| appeared, with a study of frequencies of collective 
oscillations. Their findings for the lowest breathing model are in agreement with ours. In addition frequencies of higher 
modes are calculated. 

The work was partially supported by (Spain) Grant No. FIS2005-04181, Generalitat de Catalunya Grant No. 
2005SGR-00779 and RFBR. G.E.A. acknowledges post doctoral fellowship by MEG (Spain). 



APPENDIX A: EQUATION OF STATE IN THE HIGH-DENSITY REGIME 



In the high density regime, nr^ ^ 1, potential energy dominates and properties of the system can be compared 
to the ones of a classical crystal with the lattice spacing a defined by the density a — . Several terms of the 
expansion of the equation of state can be found using harmonic approximation (for 2D dipolar crystal this has been 
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FIG. 2: Dispersion curve a;(fe) in the first Brillouin zone — ^ < fc < |- in the high density hmit nro ^ 1. Sohd hne: harmonic 
approximation, Eq. I|A2|I . dashed Une: phonons, Eq. (|A3|I . 



done in Ref. [611 )• The leading term is given by potential energy of a classical lattice 



1 f) 

Ez-W-C(3)^Kof (Al) 



3 = 



Particles move close to lattice sites Uj ~ Zj — ja <C a. The potential energy can be expanded up to quadratic terms 
in Uj and classical equation of motion can be solved looking for the wave solution Uj — roe^^'^^°'^^^\ For the dipolar 
potential this leads to the following expression for the frequency to: 

where ^(z) is the Riemann zeta function and Lin{z) is a polylogarithm fmiction. The dispersion relation (|A2p is 
shown in Fig. [51 

The low-momenta behavior |fc| — > can obtained by expanding exponents in the sums in lA2l and noticing that zero- 
and first- order terms get canceled 




uj{k) 



^ W g 1^ - 2y3C(3)^(-o)^/^|fc| (A3) 



Small |fc| behavior correspond to phonons. Same result can be recovered by calculating the compressibility in the 
system mc^ = ndfi/dn, using the leading term in the energy (jAip to calculate the chemical potential /i — dE^^^ / dN . 
Then the phononic spectra calculated as uo{k) — c\k\ exactly coincides with expression (|A3p . In Fig. [2] the phonon 
excitation spectrum is plotted against the solution (jA2[) . The latter solution has a characteristic rounding close to 
boundary of the Brillouin zone. This feature, of course, is missed in phononic description which becomes inapplicable 
for such large values of | fc | . 

Contribution i?^^' to the lattice energy (jAl[) is obtained by summation of the energy of the zero-point motion of 
atoms. Thus, one has to integrate the dispersion (|A2p over the first Brillouin zone (BZ): 



f huj{k) dk ^ , ,5/2 ,,,, 

- ' - C—^{nror^ (A4) 



N J 2 Vbz Mr'o' 

BZ 

The obtained contribution is positive, as it adds positive kinetic energy and describes displacement of atoms from the 
minimum of potential energy. Correction (jA4|) scales with density as (nro)^/^ compared to the (nrp)^ dependence of 
the dominant term. The coefficient of proportionality C = 2.26... was obtained by numerical integration of (|A4[) with 
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the dispersion relation (IA2p . 
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